Interfacing with Hamiltonian Dynamics 
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Simple constructions and protocols are demonstrated to allow the implementation of universal 
quantum computation on an arbitrarily large quantum system by controlling a fixed number of 
spins, vastly reducing the engineering requirements in comparison to a direct implementation of the 
circuit model. Analytic pulse sequences are derived for deterministic gate operation and system 
cooling which only require pulse strengths that are comparable to the intrinsic Hamiltonian of the 
system. 



I. INTRODUCTION 

The circuit model for quantum computation has been 
instrumental in the development of a plethora of exper- 
imental proposals for the implementation of quantum 
computation. One merely has to demonstrate how to im- 
plement an arbitrary single-qubit rotation on any qubit, 
and an entangling two-qubit gate between any pair of 
qubits. However, the circuit model should not be con- 
sidered a natural paradigm for the control of arbitrary 
quantum systems, merely a useful tool for proving the 
universality of the system in question. Indeed, it is often 
impossible, or extremely undesirable, to implement such 
a scheme directly. For example, single spin manipula- 
tion is extremely challenging in systems such as optical 
lattices. In other schemes, extremely high fidelity single- 
qubit operations are possible, but the number of control 
elements required for the individual addressing of each 
and every qubit in a register is prohibitive. 

This prompts the question of how little control is re- 
quired, which has been answered by proving that Hamil- 
tonian evolution, even for translationally invariant and 
rotationally invariant Hamiltonians, can implement ar- 
bitrary quantum computations [l|, l3) ISl] ■ However, these 
constructions require spins which are larger than qubits, 
or interactions with some spatial extent beyond nearest 
neighbor. Moreover, one must be able to impose suit- 
able initial conditions (a product state), which could in 
themselves require the local control that may not be 
available. Finally, error correction and fault-tolerance 
have not been built into these schemes, although there 
are some closely related schemes based around quantum 
cellular automata [j| which can be shown to be fault- 
tolerant [5|. Reintroducing a small amount of control 
into such Hamiltonian systems can markedly reduce the 
overheads involved in the constructions. One such ex- 
ample is global control, where several different Hamilto- 
nian evolutions can be pulsed, and their sequencing is 
the added control that allows the reduction to nearest- 
neighbor qubit systems in a fault-tolerant setting (see 
[3, Q and references therein) . 

In this paper, we will study another possible reduction. 



where control of a fixed number of spins is retained, and 
show how, through manipulation of these, quantum com- 
putation can be achieved in the entire quantum system. 
This is a plausible control paradigm for early experiments 
in quantum computation. Such a scheme was referred to 
as a Universal Quantum Interface (UQI) by Lloyd and co- 
authors [7[. These and others have subsequently proved 
that almost all quantum systems can be controlled in this 
way, including nearest-neighbor translationally invariant 
qubit systems Q- However, once the proof of existence 
is completed, these schemes fall back on numerical tech- 
niques for determining the control sequence for a system 
of specific size, and typically reveal little about the so- 
lution for systems of different size, or the efficiency of 
determining these sequences. Instead, we concentrate on 
slightly more complex systems where we are able to give 
concrete protocols for arbitrary system size. We are also 
able to discuss how to handle experimental shortcomings 
such as finite implementation times for the control se- 
quences. This presents a number of improvements over 
[l| , where the construction technique was first introduced 
as a tool for designing more complex Hamiltonians. For 
example, our most basic scheme still involves a nearest- 
neighbor translationally invariant Hamiltonian acting on 
4-dimensional spins, and control is over 3 of the spins. 
However, we now only need single-spin control of these 
spins, rather than needing to make them interact. We 
also discuss how these chains can be deterministically 
cooled from their initial, unknown, state, to the state 
required for the computation. 



II. QUANTUM STATE TRANSFER 

A. Basics 

We start our study by reviewing the basic components 
that we will use in Sec. IIIII to formulate the UQI sys- 
tems. We shall quote several results, with little justifica- 
tion. Their relevance will become apparent as we start 
to construct UQI schemes. The principal element here is 
the study of quantum state transfer. This was initiated 



by Bose Q, who took a simple nearest-neighbor Hamil- 
tonian on a ID chain of TV qubits 



N-l 



HsT — 2(^1 + ^n) + 2 2^ ^i^i+1 + YiYi+i, 



i=l 



assumed the chain was initiaUy prepared in the state 
|'0) |0) , and demanded how well the unknown single- 

qubit state \ip) can be transferred to the other end of the 
chain. This was facilitated by the observation that 



N 



HST,Y.^^ 



= 0, 



so the Hilbert space of the Hamiltonian decomposes into 
subspaces of fixed excitation number (number of |l)s)^. 
The state |0) is an eigenstate, and thus the task is 
to transfer from |1) |0)®^"^ to |0)®^"^ |1). In fact, for 
the purposes of designing UQIs, it is possible to entirely 
neglect the transfer of the |0) component of {ip), and just 
consider the transfer of |1). Developing the theory of 
state transfer will subsequently allow a mapping from 
the UQI Hamiltonians into a state transfer system Hst 
such that |1) |0) ~ can be considered the start of the 
computation, and |0) ^ |1) is the end of the compu- 
tation. Such a transformation can be recognized by the 
way the Hamiltonian acts on basis states, 



Hst \n) 



1> 



1) 



where, in this case, the basis states are defined to be 
\n) = |0) |1) |0) ~". In this situation, one can 

prove several basic facts, including 

Lemma 1. State transfer can be achieved with a proba- 
bility 0{N-^/^) in a time 0{N) 0y. 

Lemma 2. The arrival of {ip) = |1) can he detected by 
measuring spin N. A strategy of measure and repeat until 
the state has arrived gives a failure probability of less than 
e in a time 0{N^^^ \og{e)) [id]. 

Lemma 3. The Hamiltonian Hst — Z\ — Z^ can be 

mapped to a system of free (non-interacting) fermions. 

Given that the fermions do not interact, one can read- 
ily translate results on state transfer from the single ex- 
citation subspace to multiple excitation subspaces, pro- 
vided the phase properties of fermion exchange are cor- 
rectly handled [li[. A useful tool for this is the wedge 
product [13, IS- F'oi' example |1) A |2) denotes the pres- 
ence of excitations on spins 1 and 2. The critical property 
of this product is that |1) A |1) = 0. 



^ Note that Hst is not the Hamiltonian considered by Bose. How- 
ever, in the single excitation subspace, the two are identical. 



B. Wavepackets 

There are a number of strategies that one could utilize 
to improve the arrival of a single excitation. Haselgrove 
|14| gave a numerical technique, applicable to very gen- 
eral networks, for finding how Alice should distribute the 
excitation across the first Na spins such that it is opti- 
mally transmitted to Bob, who controls the last Nb spins 
of the chain, the idea being to work in the regime where 
Na -\- Nb ^ N. However, consideration of the special 
case of a linear chain, enables the derivation of rigorous 
lower bounds, utilizing the following fact: 

Lemma 4. In an infinite chain, away from the ends of 
the chain, the optimal encoding of a single excitation is 
in a Gaussian wavepacket of spread A ~ N^^'^ , where 
N is the distance of transmission [15]. 

The wavepacket that has the maximum group veloc- 
ity for the system, and is minimally dispersive, can be 
written using the previously defined position basis as 



IV'g) = 



1 



VAV^„f 



E 



^-{n-xoy/(2A'')-i7Tn/2 | x 



(1) 



which is centered on spin xq at time t — 0, and its center 
is located around Xq -\- 2t at any later time t. The speed 
is dimensionless due to the choice of units of, h — 1, and 
the suppression of the energy scale of Hst. 

In addition to the deterministic arrival and low disper- 
sion of the wavepacket, a further advantage is that the 
wavepacket is well localized. This means that multiple 
excitations can be sent at a rate of 0(N/A), without 
degradation of the signal [12, [l3| ■ 



C. Time Control and Wavepacket Injection 

While the optimal wavepacket for transmission 
through the central region of a long chain is known, this 
is of little benefit since we must find the encoding at the 
end of the chains, as well as needing to find how to inject 
the wavepacket onto the chain, and remove it at the end. 
To this end, consider a second situation introduced in 
[14J , where time- varying control at the ends of the chains 
is available, 

2Htv = no{t)iXoYi + XoYi) + r!Ar(i)(Xjvlw+i + XnYn+i) 



7 ^ XiXi^i + YiYi+i. 



This time-varying control can be used to simulate an in- 
finite spin chain i.e. to make the N spins of the chain 
behave as if they were in the middle of a much longer 
chain, introduce a Gaussian wavepacket on one end, and 
remove it from the other end at a time 0{N) later. The 
state of the chain to be simulated, as a function of time. 
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FIG. 1: Plot, as a function of time, for (a) Fidelity of arrival and departure of states |2)® and [0)®^ respectively for / = 3 
and A — 10, giving an arrival fidelity of 99.45%, and (b) coupling strengths Qq and fl;v. The final fidelity can be increased by 
increasing A. 
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Similarly, the state of the actual chain is 

<p7i{t)\n) , 



N+l 

E 



where the aim is to ensure that V'n(i) — 4'n{t) for alH > 
and n = 1 ... TV. At i = 0, the packet starts in a nega- 
tive position, xq = —/A, where / encapsulates an error 
parameter determining how many standard deviations of 
the wavepacket we wish to use, ultimately influencing the 
error of the simulation. In the time t = /A + N/2 the 
wavepacket moves along the infinite chain with a ballis- 
tic motion, ending up at a position N + /A, outside the 
region of interest. Using Eqn. ([1]), 



^-(ri-3:o-2i) /(2A^)-i7rn/2 



In order for the actual chain to simulate the infinite chain, 
it suffices to ensure that V'n(O) = 0n(O) and -^ — -^ 



dt 



for all times. This is convenient because the Hamiltoni- 
ans are identical in the two cases, except at the ends of 
the chain. Thus, we have 



from which 



i—7- = f^o0o + 02 = -00 + "02, 
at 



O - ^0 

\Iq — -— 
00 



is readily derived. It is worth noting that |0oP = 
Sn=-oo IV'nP, which could be used to determine |Slo|. 
However, an alternative derivation^ (since the system is 



^ In the limit of A -^ oo, the two derivations give the same 

pulse. In fact, for small A, this derivation, by setting \4>o\'^ = 
1 

In=-oo l'/'(^)P'^^ii seems to give better results numerically, but 
the pulses extend over longer time 



not continuous), which allows the modulus sign to be 
removed, derives from 



t-^ = noit)Mt) ^ ¥^ = Mt)Mt), 



which ultimately reveals that 

no{t) = ^= 



VAV^Wl-5e-i/4A^ 1 + Erf 



./ 2t+ao-^ 



where Erf(z) is the error function. At i = 0, flo{t) « 0, 
and for all times greater than 0(A) (i.e. once the bulk 
of the wavepacket is on the chain), this is true again. In 
the meantime, flo{t) is a smooth function that is never 
greater than 1, as can be seen from the numerical exam- 
ple depicted in Fig. [TJ The total time for the sequence is 
0{N + N^^'^loge), where the second term incorporates 
the influence of increasing the accuracy by increasing /. 
Similarly, one can derive that 

flNit) = ^= 



VaT^I 1 - ie-V4AM 1 _ Erf 



A 



The boundary condition that is imposed on 0jv+i(O i^ 
0Ar+i(oo) = (— i)^+^, which is due to a phase issue that 
occurs during a state transfer protocol (see, for example, 

M)- 

Although concentrating on the specific task of trans- 
ferring an excitation, this analytic solution is applicable 
to the transfer of an unknown single qubit state lip) ini- 
tially stored on spin onto spin A^ -|- 1 simply due to 
linearity - the |1) state transfers, and the state contain- 
ing no excitations is an eigenstate. In the remainder of 
this subsection, we will show that the protocol is very ro- 
bust against control errors and, unlike other state trans- 
fer schemes, the available control allows the preparation 
of the initial state of the system rather than having to 
assume an additional cooling mechanism. 



1. Robustness 

The described protocol usefully indicates how a num- 
ber of very relevant physical effects can be incorporated 
into the scheme, such as finite laser power. There are 
two further experimental concerns, decoherence and im- 
perfect control operations. We shall start by numeri- 
cally testing the effects of imperfect control operations. 
Take, as an example, the case of A^ = 100, A = 30 and 
/ = 3. To simulate the time evolution, we split the se- 
quence into time steps of 0.1 (compared to the speed of 
the wavepacket, which is 2), and assume that the values 
of r^o.Af (i) are constant over that timescale. For the ideal 
pulse sequence, the success probability is 99.94%. This 
probability tends to 100% as / -^ oo and A/N^^^ -^ cx). 
Even after failure, all is not lost because we can con- 
tinue to perform a heralded state transfer protocol until 
successful. 

First, we consider timing errors i.e. Cli{t) = fli{t — ti). 
Evidently, to determines when the protocol starts, and so 
the only relevant parameter is the relative time tjy — to- 
A plot of calculated values is given in Fig. [3] 

The next imperfection that we consider is the strength 
of the laser pulse. For example, with a systematic offset 
of 5%, i.e. Cti{t) = 1.05fli{t), the arrival probability is 
still 99.2%. Finally, we consider fluctuations of the con- 
trol fields, by a random amount at each time step, where 
we considered two different cases; where the maximum 
fluctuation was up to 10% of the ideal field at that time 
(arrival probability 0.9988(1)), and by up to 0.01 (com- 
pared to the maximum field strength of approximately 
0.3), independent of the value of the intended field at 
that time (arrival probability 0.9986(2)). An example of 
one instance of this evolution is given in Fig. [21 These 
results suggest a large degree of tolerance to control er- 
rors. 



2. System Initialization 

In addition to using the temporal control for injecting 
wavepackets into the system, it can be utilized to provide 
an effective cooling from an unknown initial state to the 
state |0) . We shall assume that we have the ability to 
measure spins and A^ -I- 1, and to reset them to the |0) 
state. As before, the time-varying controls at the end of 
the chain can be used to simulate a larger chain, designed 
such that, with high probability, after a fixed period of 
time, excitations leave the part of the chain that they're 
initially localized on. We start by designing a protocol 
under the assumption that there is a single excitation 
somewhere on the chain. Working sequentially, starting 
with spin 1, pulse sequences can be designed so that if an 
excitation is there, it moves off onto spin or spin A^-|- 1, 
and can be reset. When spin is measured, it reveals 
whether the excitation was removed or not. If not, the 
protocol continues. If the excitation had originally been 
on spin 2, it is no longer localized there. However, the 



state that it would be in after the first pulse sequence 
can be found via an efficient classical computation. So, 
a pulse sequence can be designed to remove it. This 
continues for all N spins, after which the single excitation 
is guaranteed to have been removed. 

In the case where multiple excitations are present, ex- 
actly the same protocol serves to remove all excitations, 
except that we do not stop the protocol once a single exci- 
tation has been removed. For example, let's assume that 
there were originally 3 excitations, one on each of spins 1, 
2 and 4, which is denoted by |1)A|2)A|4). The first 4 pulse 
sequences that we go through can be described by Ui to 
U4 respectively. After the first evolution, Ui\l) = |0), 
which is detected and removed, leaving Ui |2) Af7i |4) (es- 
sential to this is the fact that (0| C/i |2) = (0| t/i |4) = 0). 
After the second evolution, we have 

U2U1 12> A U2U1 |4) = |0) A U2U1 |4) , 

Again, since (0| U2U1 |4) — 0, this is just equal to 
U2U1 |4). After the third pulse is used, no excita- 
tion is detected on spin 0, but on the fourth we have 
U4U3U2U1 |4) = |0), and all the excitations are success- 
fully removed. 

All that remains is to specify the chain to simulate 
with fli{t). One str ateg y is to use the perfect mirroring 
chain introduced in [iTj, where it was proved that under 
the action of the Hamiltonian 



HpST = 5 X! 

n=l 



^/n{M-n){XnXn+l + YnYn+i) 



for a time 7r/2, an excitation initially on spin n arrives 
perfectly on spin M + 1 — n. For a chain with a large 
number of qubits M, the central iV spins are approxi- 
mately uniformly coupled, with a value of M/2. Thus, 
a good approximation for simulation of this system can 
be achieved, after rescaling to 2Hpst/M. By selecting 
the mirroring point to be just offset from the center of 
the chain, excitations get moved off the ends. The exact 
solutions for the functions ipmit) = {m\ e~*-^^sTt j^^^ ^^^^ 
be found in 11] and, as such, we are able to take the limit 
of A/ — > cx) to achieve an arbitrarily accurate simulation. 
To prove this, consider the perturbation 

(M+N)/2 

6H^ Y, {M/2-y/n{M - n)){\n) {n + l\+\n + 1) (n|) 

n={M-N)/2 

acting on HpsT- By ensuring that 5H retains the nur- 
ror symmetry of the original Hamiltonian, the analysis 
is simplified because it is only the eigenvalues that affect 
the fidelity of state transfer pjQJ. Provided the matrix 
elements of 6H are much smaller than 1, i.e. M ^ N'^, 
perturbation theory allows us to calculate how the eigen- 
values are shifted, using (A„| SH |A„), where |A„) are the 
unperturbed eigenvectors. The error can be bounded by 
using the eigenvector which has the largest elements in 
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FIG. 2: Plot, as a function of time, for (a) Fidelity of arrival and departure of states |2) and |0) respectively for / = 3 

and A = 10, giving arrival fidelity of 99.37%, and (b) coupling strengths CIq and CIn, which are randomly varied from their 
ideal forms (Fig. [l| within the range ±10%. To be compared with Fig. [l] 



the middle of the chain, 



|Ao) = 




Taking the large N limit and using Stirling's approxima- 
tion, we find that 



(Ao|<5i/|Ao> 



7V3 

M572- 



Provided this is small, all higher order terms from per- 
turbation theory are even smaller and can be neglected. 
Small errors of size 5 in the eigenvalues contribute to 
a reduction in arrival fidelity of the state mirroring of 
0((5^) [IJl, and thus we are guaranteed that in the large 
M limit the error tends to 0. 



D. A Second Time-Controlled Scenario 

There is a second instance that we will be interested in, 
where the scheme is no longer a perfect chain, but the net- 
work depicted in Fig. Ufa). Again, the aim is to achieve 
perfect excitation transfer between the shaded spins, or 
to place a lower bound on how well this can be achieved. 
We would aim to do this, as before, by simulating a much 
longer chain. However, when trying to transmit a Gaus- 
sian wavepacket through the system in Fig.[51[b), part of 
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FIG. 3: Arrival fidelity as a function of timing error, ijv 
for a chain of 100 spins, A = 30, f ~ S. 
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the packet gets reflected at the T-junction, and trans- 
mission is no longer perfect. A heuristic explanation for 
how to get around this problem is that we know that 
wavepackets can be created that travel along the chain 
with speed 2. When they hit the junction, they are par- 
tially reflected. If the side that Alice controls is short 
(i. e. Na -^ oo), and the side that Bob controls long, 
then the first reflection will rebound off Alice's end, and 
impinge once again on the T-junction, part of which will 
get transmitted. Everything that gets transmitted can be 
received by Bob, provided his region of the chain is long 
enough that the wavepacket does not reflect and move 
out again. If the transmission probability is T, then after 
TO such cycles (requiring a time tuNa) some proportion 

which is roughly X]fc=o -^(^ ~ ^)* ^^^ passed beyond the 
imperfection, and can be gathered by Bob after a further 
time N/2, provided Ng > 2{m~ 1)Na- For a probability 
of failure e, we would only require to, ^ loge provided T 
is a constant'^. 

The first step towards proving this is to show that for 




FIG. 4: (a) Time-controlled spin chain, used to simulate the 
fixed spin chain depicted in (b) , where . . . indicates the pres- 
ence of a chain of Na (left) or Nb (right) qubits, either of 
which may be taken to oo. 



After calculating the variation of T as a function of the energy, 
this can be derived more rigorously, taking into account the fact 
that the reflected wavepacket is no longer Gaussian, and therefore 
its transmission probability is different. 



a half-infinite chain, when a wavepacket is incident on 
the end, it is reflected as a wavepacket that continues 
to evolve as if the boundary was not present. To see 
this, consider an infinite chain where the initial state is 
an equal superposition of two wavepackets (taken, for 
simplicity of exposition, to be Gaussian) positioned at xo 
and —xq, traveling towards each other with equal speeds 



V, 



i\^Gixo,-v)) - \^:Gi~xo,v)))/V2. 



By linearity, it is clear how these packets evolve; they 
just move smoothly past each other and, at long times, 
are two Gaussian wavepackets moving away from each 
other at equal speeds, 

i\ipG{xa - vt, ~v)) - \ipG{~xa + vt, v)))/\f2. 

The Hamiltonian of the infinite system in the single exci- 
tation subspace commutes with the symmetry operator 



^ = Ei 



|0> (0| 



and, as such, the Hilbert space can be decomposed into 
two subspaces, symmetric and antisymmetric. The de- 
scribed wavepacket is antisymmetric, and thus its evolu- 
tion is described entirely within that subspace, which is 
precisely that of a Gaussian wavepacket on a half-infinite 
chain impinging on the boundary. Thus, we learn that 
wavepackets reflect perfectly'', although they acquire a 
phase of tt. It still remains to prove the expected reflec- 
tion off the T-junction. Having proved what happens at 
the boundaries, it is now sufficient to take the size of the 
simulated systems N^ , Ng — > oo , and show what hap- 
pens when a wavepacket is incident on the imperfection. 
This proof will closely follow the work of Childs, [18| . 
The first step is to find the eigenvalues and eigenvectors 
for the network in Fig.UJ^b). Ignoring the T-junction, the 
unnormalized eigenvectors would take the form 



E 



e-^'^'^ \x) 



X— — oo 



x— — c« 



x=l 



and solve the set of linear equations such that the en- 
ergy is still 2cosfc. This reveals that the transmission 



* By which we mean that when the wavepacket is a long distance 
from the boundary, it is Gaussian 



probability for a wave with momentum k is 



4sin2(2fc)(Ji2 



J|- 4 cos2fc)2 



4sin^(2A:)(j2 + J^ - 4cos2 k) 



(J|-4cos2fc)2 



By selecting J2 = Acos^ k, perfect transmission can be 
realized for a specific value oi k. To use the minimally 
dispersive wavepacket, with maximum transmission rate, 
k = tt/2 is the ideal choice, so J2 = is selected (i.e. rip = 
Ji is constant). Expanding for small 6k about this point, 
we find that 

4 



l^7r/2+(5fel 



1 



74 
•^1 



5k^ + 0{5k^). 



For a Gaussian wavepacket of width A, the width in k- 
space is roughly 1/A, and thus the transmission proba- 
bifity -- 1 - OiN-"^/^), so in the fimit of large iV, the 
success probability tends to 1. In the instance where 
only the right-hand side of Fig. HKb) is semi-infinite, and 
the left-hand side only extends to ~Na, any component 
that is reflected off the T-junction is perfectly reflected 
off the flnite end, and is subsequently incident on the T- 
junction again. By waiting long enough to collect m such 
repetitions, the success probability is enhanced to 



1 



j4A"+i(m + l)! 



O 



1 



/\^rn+3 



with energies 2cosfc. However, the T-junction causes 
partial reflection of the wave. Thus, we write the new 
eigenvectors incorporating the amplitudes for transmis- 
sion and reflection, 

—1 00 

J2 (e-*"+i?fee''=")|:r)+^rfce-^'=^|x)-f ^ a,\i) , 

ie{0,Q,/3,7} 



The choice of m is related to the length of time spent 
collecting the wavepacket, m ^ {t/2 — N) / {2N a) ■ Thus, 
the total scaling is also 0{N -\-N'^''^ loge) for a maximum 
error of e, and the entire wavepacket still has a width 
(9(7V^/'^), so the scaling of the rate of transmission is the 
same as the previous scheme. 



III. DESIGNING A UQI 

Having stated and derived all the basic tools that are 
required for the construction of a UQI within the frame- 
work of state transfer, we are in a position to start de- 
signing Hamiltonians whose dynamics can lead to quan- 
tum computation when aided by dynamical control over a 
fixed number of spins, by proving mappings to the state 
transfer results. In order to prove universal quantum 
computation, it is sufficient to show how to perform any 
single-qubit rotation on any qubit, and how to perform 
just one entangling gate between an arbitrary pair of 
qubits. Given that we retain control over at least two 
qubits, in principle it suffices to show how to place any 
pair of qubits on those two control qubits, and the arbi- 
trary control over them allows the implementation of the 
desired gates. So, all we need is a simple permutation op- 
eration. This motivated the original UQI construction in 
[l|, acting on a chain of N 4-dimensional spins, whose lo- 
cal Hilbert space is decomposed into two qubits, denoted 
a and 6, 



JV-l 



Hi 



2 z_^ 



{XX + YY)l 



n,n+l 



'SWAP 



n,n+l ■ 



To see how this works, start by observing that 



N 



0, 



and adopt a basis in the single excitation subspace of the 
a qubits of |n), n = 1 . . . A'^, denoting the presence of the 
single excitation on spin n. The action of Hi on these 
basis states can be written as 

ffi In)'^ \i,S = \n ~ ir IV-n-i)' + I" + 1>" IV-n+i)' , 

where \ipn+i) ~ SWAP„_„+i \ipn), and thus this Hamil- 
tonian maps into the state transfer system Hst ^ Zi — 
Zn. With control over spins 1 and N, Z\ and Zat 
can be incorporated. Thus, by invoking Lemma [1] if 
the system were to start in |1)° j-i/;) , then after a time 
0{N), with probability 0(7V-V3), j|. ^0^1^ be in state 

|A^)"(SWAP^,_l,^,...SWAPl,2|V'))^ which is a cyclic 
permutation of the h qubits. Multiple repetitions of this 
protocol enable any arbitrary qubit to be placed on spin 
1. Moreover, by measuring to see if an excitation is on 
I A^) , the arrival is heralded, and the protocol of Lemma 
|2] gives arrival with failure probability less than e in a 
timeO(iV5/3log(e)). 

Similarly, we could start an excitation on |2) j-f/;) 

and transfer it to liV)" (SWAPtv-i.at • . . SWAP2,3 |V'»\ 
which also implements a cycling operation, but only on 
qubits 2 to N . Thus, any pair of qubits can be placed 
on the first two spins by controlling only the first 2 spins 
and spin N . In order to implement a one or two-qubit 
gate requires a time 0{N^I'^\o^(e)). 

There are now a number of potential improvements 
that we would like to introduce. Firstly, we might like 
to restrict our control of spins to being strictly local, i.e. 
no two-qubit interaction. We do this by recalling a re- 
sult of [IQI , whereby iV -I- 1 repetitions of global applica- 
tions of controlled-phase {CP^ gates and Hadamard gates 
(iJ) gives the mirroring operation between the N qubits 
(i.e. the state of spin n is mapped to that of A^ -t- 1 — n). 
By altering Hi , replacing SWAP with the new operation 
U = {ii®l)- CP, we get 

Af-l 

H,= Y,\Ql){lQ\l^^,®U + h.c. 



Note that because U ^ U\ it has become necessary to 
split the |(XX+Fy) in order to maintain the Hermitian 
nature of the Hamiltonian. We interpret this as being 
that when the excitation steps to the right, it implements 
U on the 6 register, and when it steps to the left, it im- 
plements C/t ^ which means that the position of the exci- 
tation alone is sufficient to reveal what sequence of gates 
has been implemented. Now when we perform a heralded 
state transfer protocol from 1 to N , the implemented uni- 
tary is [/ = {Hi®. . .®Hn-i)[CPn-i,n ■ --GPia), which 
is the aforementioned global application of controUed- 
phases followed by Hadamards, except that a Hadamard 



needs to be applied on the h qubit of spin N . Thus, N -\-\ 
repetitions of the protocol implements the mirroring op- 
eration on spins 1 to A^. Similarly, performing transfer 
from 2 to A'^ implements mirroring on spins 2 to N . The 
combination of the two is the previous cyclic permuta- 
tion; thus we recover the ability to place any single qubit 
on spin 1, although it now means that a one-qubit gate 
requires time 0(A^^^/"^ log(£)). To achieve a two-qubit 
gate, we run the state transfer protocol from 1 to A^, ap- 
ply a Hadamard gate on qubit l'', and then run the state 
transfer protocol from A^ to 1, which implements the in- 
verse unitary. Thus, U^ HU is a two-qubit gate acting on 
spins 1 and 2 which is locally equivalent to controUed- 
NOT. So, we have an entangling gate between spins 1 
and 2, and cycle. However, the controUed-NOT can be 
composed to give SWAP, so that as the qubits are cycled, 
SWAP can cause a single qubit to stay on spin 1 while the 
others cycle. Therefore, this can be used to implement 
the two-qubit gate between any pair of qubits. 



A. Realistic Pulses 

So far, we have assumed that we can place excitations 
on spin 1, remove them from spin N etc., arbitrarily 
quickly, at least compared to the inverse of the Hamil- 
tonian interaction strength. This is not a very realistic 
condition, and one that can be remedied. One of the 
problems with a slow laser pulse which converts |0) to 
|1) on a particular spin is that when it has been partially 
rotated, the Hamiltonian is already moving that part of 
the excitation onto other spins. As a consequence, it is 
possible that the pulse would introduce multiple excita- 
tions. One trick to avoid this involves introducing a third 
level, 1 2), to system a such that when an excitation hops 
to the right due to the Hamiltonian, it leaves behind it 
a 1 2) rather than a |0). Thus, a laser J7o(|0) (1| + |1) (0|) 
cannot create excitations on the |2) component. Such a 
Hamiltonian takes the form 



JV-l 
n=l 



(101 



n,n+l 



c^+ioi>(i2i:,„+i' 



h.c 



When enhanced by the laser fields 

Hl = f^o(lO) (1| + |1> (0|)? + nNi\2) (1| + |1> (2|)^, 

we can now prove the mapping to the state transfer 
scheme of Htvj so |0) is coupled to |1) |0) ~ with 
strength flo{t), and the central states are of the form 
|2)®""^|1)|0)®^"". Finally, the f^^v pulse leaves the 
system in the state |2) , from which it must be re- 
set before the next U can be applied. This is achieved by 
applying a laser on spin 1 to convert the |2) into a |1), 
and on spin A^ to convert a |1) into a |0). The second 
term in H3 then takes care of the rest. 

For the part of the sequence where we start the exci- 
tation on spin 2, the scheme maps into that of Sec. IIIDI 



One can see this because we apply a laser pulse to the 
state |0) , which creates the state |0) |1) |0) ~ , and 
the two are coupled by the laser strength fli{t). The 
Hamiltonian further couples this state to two states, 
|1)|2)|0)®^"' and |0) |2) |1) |0)®^-^ the latter then 
couples further, forming a chain of states. With a third 
laser on spin 1, couphng the states |1) and |2), the 
state |1) |2) |0)®^~' is coupled to |2) |2) |0)®^-' with the 
laser strength flo{t). Consequently, the state transfer 
sequence requires 0{N + N^^^loge), and thus a single 
gate, one- or two-qubit, can be achieved within a time 
0{N'^/^ + N^/^loge), which is more efficient than the 
previous 'fast' pulses, and derives all the benefits of ro- 
bustness demonstrated in Sec. IIIC II 



B. Cooling 

Naturally, it is desirable to be able to prepare the ini- 
tial state of the system using the control available. With 
the described mapping between the UQI system Hl + H3 
and a state transfer chain Htv, one would hope that the 
first step in preparing the initial state of the system would 
arise from application of the results in Sec. IIIC 21 The 
mapping would allow us to remove any of the |1) states 
from the a system. Unfortunately, this is not the case 
- we have only proved the mapping in the first excita- 
tion subspace. In order for the results to be applicable in 
higher excitation subspaces, one of two conditions needs 
to be fulfilled. Either the rotations U need to commute 
with each other (i.e. [Ui^is,Ui^i+i] — 0) or the state of 
the b system would have to be in the maximally mixed 
state p = 1/2^. While this second condition may often 
be a good assumption, we note that a more generally ap- 
plicable technique was described in [8| where one simply 
measures one of the spins, and sets it to a fixed state (say 
|0)), and repeating. Note that unlike the rest of [9], it is 
not necessary to use an ancillary system. 

Once all the |1) excitations have been removed from 
the a subsystem, the entire system is in an eigenstate of 
the Hamiltonian. However, we still have to prepare the a 
system as |0) rather than an unknown superposition 
of |0) and |2) states. Consider the setting where all spins 
are |0) except for one, which is in the |2) state (this ar- 
gument trivially generalizes to any arbitrary distribution 
of |2)s). Now, each time wc apply the standard state 
transfer protocol from spin 1 to spin N , the |2) shifts one 
spin towards spin 1 (and flips to a |0), where all other 
states have become |2)). Thus, after no more than N 
repetitions, the unwanted state has reached the control 
area and can be corrected. Once the a spins are cor- 
rectly initialized, we can implement the cycle operation 
on the b qubits, setting each of them in turn to |0) as 
they pass through a control spin, ready to proceed with 
the computation. 



Fault Tolerance 



The architecture that we have described so far does not 
support fault-tolerant quantum computation. With only 
a fixed number of access points to the chain, the degree of 
parallelism is insufficient, which means that errors build 
up more quickly than they can be corrected. One can 
avoid this problem by allowing control of a regular array 
of spins [3]. However, this is not the end of the problems. 
We will not dwell on these here, merely describe them, 
and possible solutions, for the sake of completeness, while 
noting that much of the structure of the scheme (limited 
parallelism, the need to correct both a classical and quan- 
tum region of the device, when the only way to interact 
with the quantum region is through the classical region) 
is very similar to global control schemes, in which fault 
tolerance is possible Q. Moreover, the present setting 
has the advantage that we can perform different actions 
(e.g. error correction based on syndrome measurements) 
at each point of interaction with the chain, rather than 
needing the same interaction. As such, one might expect 
a vastly improved fault-tolerant threshold. 



During the computation, no matter how well we pro- 
tect the system, some noise will always enter. Typically, 
one assumes that these errors are single-spin rotations. 
Of course, if a single b qubit is affected by noise, then 
standard error correction techniques are sufficient to cor- 
rect the error. However, if an error were to flip, say, 
between the |0) and |2) states of an a spin, this means 
that every time a global unitary is applied through the 
state transfer protocol, it is faulty on a specific pair of 
qubits. As with the cooling mechanism, this pair changes 
with each repetition, so that by the time the error has 
propagated to the control area, many of the computa- 
tional qubits in the block will have been affected. The 
effects are somehow related, and it may be possible to 
design error correcting codes to compensate for this. In 
particular, our error detection sequence only needs to de- 
termine where the flip first occurred, which is very sim- 
ilar to detecting the edge of a pattern of spins, such as 
the boundary term in |11 . . . 100 . . .0). This can be re- 
lated via a local unitary transformation into detection of 
a single excitation [20|, which is precisely the task that 
standard error correcting codes are designed for. Alter- 
natively, given that the error only persists for a finite 
time, then the affected region of the b qubits is no larger 
than the separation between two control regions, and er- 
ror correcting codes that are capable of correcting that 
many errors can be used (or the error correcting codes 
can be used on non-contiguous blocks of qubits). A sim- 
ilar effect will result if the initial Hamiltonian is subject 
to local perturbations, or if multiple |1) states were to 
appear in the a region. Thus, it should be possible to 
handle all such issues, although it is not our priority in 
the present paper. 



IV. SUMMARY AND OUTLOOK 

We have presented a constructive technique that shows 
how simple one-dimensional Hamiltonians can be used for 
universal quantum computation, with control of only a 
small number of the spins. This construction is based 
on using state transfer in an auxiliary system and thus, 
when limited to a nearest-neighbor translationally invari- 
ant Hamiltonian necessarily requires a local Hilbert space 
dimension of at least 4 (2 for the auxiliary system, and 
2 for the computational qubit). While control is over 3 
spins of the system, each of the actions we require is re- 
stricted to acting on a two-dimensional subspace of an 
individual spin. Furthermore, we have shown how to in- 
corporate a number of useful aspects such as robustness 
to control errors, a deterministic cooling protocol and the 
slow nature of the control fields. 

Evidently, the intention of designing such a scheme 
is to make experimental implementation more feasible. 
However, more accurately, this work represents a limit of 
analytic constructive techniques, illustrating the abilities 
that are present in simple systems. The first experiments 
will require even simpler systems, such as a ID Heisen- 
berg chain of qubits, and, as such, one expects the con- 
trol theory solutions to be the first to be experimentally 
tested, and the results here only become relevant when 
scaling becomes a critical issue in the design of control 
pulses in these simpler systems. 

An interesting avenue for further theoretical work 
would be Hamiltonian simulation - our constructions 



have been carefully tailored to give quantum gates that 
can be composed using the circuit model. However, one 
of the major applications of a quantum computer will 
be quantum simulation, where we attempt to simulate 
a variety of different Hamiltonians in a system which is 
easily controlled. Are there UQI systems where control 
over the spins influences an effective Hamiltonian evolu- 
tion, without having to resort to a Trotter decomposition 
of the evolution? If resorting to a Trotter decomposition, 
one solution immediately suggests itself. Temporarily ne- 
glecting the 'slow rotations' component, and taking the 
a system to have 4 levels, a Hamiltonian 



^sim = EE|Oz)(*0|;„^l®C/^ 



h.c 



where Ua = e~*'^*°''^°' and o-^ G {X,Y,Z}, allows simula- 
tion of a Hamiltonian J^ni'^x-^-^ + \^^ + ^zZZ)n,n+i 
simply by selecting the correct ratio of the number of uses 
of the different excitations in the state transfer process. 
Furthermore, by applying the inverse (i.e. by implement- 
ing state transfer on the a system from spin N to 1), 
we are able to generate composite pulse sequences to in- 
crease the accuracy of the simulation [21| . 
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